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ABSTRACT 



We present an analysis of five X-ray Multi-Mirror Mission {XMM ) observa- 
tions of the anomalous X-ray pulsar (AXP) IE 2259-1-586 taken in 2004 and 2005 
during its relaxation following its 2002 outburst. We compare these data with 
those of five previous XMM observations taken in 2002 and 2003, and find the 
observed flux decay is well described by a power law of index —0.69 ± 0.03. As 
of mid-2005, the source may still have been brighter than preoutburst, and was 
certainly hotter. We find a strong correlation between hardness and flux, as seen 
in other AXPs. We discuss the implications of these results for the magnetar 
model. 

Subject headings: pulsars: individual (IE 2259-1-586) — X-rays: stars — stars: 
neutron 



1. Introduction 



It is now commonly believed that soft gamma-ray repeaters (SGRs) and anomalous 
X-ray pulsars (AXPs) are n eutron stars with ultra-strong magnetic fields, i.e. magnetars 
( jPuncan fc Thompson! Il992l ). Their common nature was conclusively demonstrated when 
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AXP IE 1048.1-5937 was observed to emit SGR-like bursts in 2001 flGavriil et all l2002h 
and IE 2259+586, in the su pernova remnant (SNR) CTB 109, wa s seen to undergo a major 
SGR-like outburst in 2002 jKaspi et al-lboosl : IWoods et al.ll2004l ). Subsequently, a variety 



variations (Gavriil & Kaspi 2004; 


Dib et al. 


2007) and slow and rapic 


pulse profile chaneres 


(Iwasawa et al. 


1992 


Kaspi et al. 


2003: Woods et al. 2004: Israel et al. 


2007: 


Dib e 


t al. 2007. 


2008 


), in addition to bursts and outbursts ( 


Gavriil et al. I2OO6: Woods et al. 


2005: 


Dib et al. 


2007 


: see Kaspi 


2007 


for a recent review). 



During IE 2259+586's 2002 outburst, the pulsed and persistent fiuxes rose suddenly by a 
factor of >20 and decayed on a timescale of months. Coincident with the X-ray brightening 
the pulsar suffered a large glitch of fractional frequency change 4 x 10~^ (iKaspi et al.l 12003 



Woods et al.ll2004l ). In the first few hours of the outburst, the pulsar's pulse profile changed 
significantly, its pulsed fraction decreased, and its spectrum hardened drama tically. Over 



80 sh ort SGR-like bursts from the pulsar were observed at the same time (IGavriil et al. 
2004r) . A near-infr ared (K^) enhancement was also observed during the epoch of the outburst 
kaspi et al.ll2003l ). 



Combining Rossi X-Ray Timing Explorer (RXTE ) observations a nd XM M observations 
of IE 2259+586 taken during and after the outburst, IWoods et al.l (120041 ) found that the 
decay of IE 2259+586's unabsorbed fiux (mostly inferred from RXTE pulsed fiuxes) after 
the outburst was well characterized by two power law components: a rapid steep decay 
visible only during the first several hours ( < 1 day) of the o utburst, and a slower decay of 
index —0.22 for the next several months. iTam et al.l (120041 ) found that the near-infrared 
enhancement at late times decayed at the same rate as the slow X-ray decay, although there 
were no IR observations during the first few hours of the outburst. 

Other AXPs have also exhibited transient behavior that could be explained by SGR- 
like outbursts. AXP XTE J1810— 197is called transient because it was only discovered in 



2003 w hen it suddenly beca n ae bri ghter by a factor of 100 (jibrahim et al.ll2004l : iGotthelf et al. 



2004[ ). iGotthelf fc HalpernI (120071 ) found that the fiux of XTE J1810-197after 2003 followed 
an exponential decay of timescale 233.5 days. Similarly, the AXP CXOU J164710. 2— 455216 
was found to have brightened by a factor of ~300 between two XMM observations taken 
5 days apart in 2006 September Jlsrael et al.ll2007h . Candidate AXP AX 1845-0258, was 



discovered in an observation made in 1993 by ASCA (IGotthelf fc Vasishtlll998l : iTorii et al. 



19981 ) . Follow-up observa ti ons in 1999 showe d that the source's fiux was smaller by a factor of 
~10 JVasisht et allboooh . IXam et all J2006h found that AX 1845-0258 remains undetected 
in C/iandra observations taken in 2003, with its fiux ~260-430 times fainter than observed 
in 1993. 
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The transient AXP phenomena summarized above are quahtatively similar to the 1998 
August 27 flare of SGR 190 0+14, in which the X-ray flux decayed with a power law of index 

0.9 JPeroci et al.ll2003h . and the flux decay of SGR 1627-41 since 1998, whi ch followed 

a power law of index ~ —0.47 and lasted for ~800 days (IKouveliotou et al.ll2003l ). However, 
thus far, the AXP outbursts have been much less energetic than most SGR outbursts. Also, 
most of the burst energy was released during the afterglows of the AXP outbursts, while for 
SGR outbursts, the X-ray afterglows have less integrated energy than the burst itself. 

With now a handful of AXP and SGR outbursts and subsequent relaxations observed, 
we can begin to look for correlations between different outburst and relaxation properties 
in the hope of constraining magnetar physics. For example, SGR outburst recoveries have 
been modeled as crustal cooling following impulsi ve heat injection, and i n principle can yield 
constraints on the nature of the crustal matter (ILyubarsky et al.l 120021 ). Alternatively, the 



AXP events have been interpreted in terms of magnetospheric twisting (IThompson et al. 



2OO2I : iBeloborodov fc Thompson! 120071 ). whose recovery depends on electrodynamics in the 
region of th e magn etosphere immediately outside the stellar surface. On the other hand, 
Giiver et al.l (120071 ) suggest that AXP recoveries can be modeled with a stationary magne- 
tosphere, with only the surface temperature changing. They argue that their model, which 
includes the stellar atmosphere, can be used to quantitatively determine the source's mag- 
netic fleld. 

In this paper we present a spectral and pulsed flux analysis of 10 XMM observations 
of AXP IE 2259+586 taken between 2002 and 2005, as the source relaxed back toward 
quiescence following its 2002 outburst. We compare the X-ray flux and spectral evolution 
of IE 2259+586 with those of other magnetars, and interpret these results in terms of the 
magnetar model. 



2. Observations 



2.1. XMM-Newton Observations 



Ten XMM ( iJansen et al.ll200ll ) observations were analyzed for this paper. The flrst flve 
observations of IE 2259+586 were taken between 2002 and 2003, just prior to and after 
the 2002 Jun e outb urst. Data from these flve observations have already been presented in 
Woods et al.l (120041 ). We re-analyzed these observations using the XMM calibrations pub- 
lished on 2007 September 4 (XMM-CCF-REL-23dJ. The later flve observations were taken 



http: / /xmm.esac.esa.int / external/xinm_sw_cal/calib / reljiotes/index.shtml 
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between 2004 and 2005. Most of these observations pointed at IE 2259 +586, with the 
European Photon Imaging Camera (EPIC) pn camera (jStriider et al.ll200ll ) in Small Win- 
dow Mode. However, three observations were obtained with XMM pointing at a portion of 
the SNR CTB 109's shell and with the pn camera in extended Full Frame Mode. Details 
about the observational mo des, pointing offsets , and exposure times are presented in Table 
[H The EPIC mos cameras (ITurner et al.ll200ll ) were operating in Full window mode with 
the medium filter in four of the first five observations, the exception being the third, and 
therefore the observed spectra are highly piled-up. The mos cameras were operated in Small 
Window Mode with the thick filter in the remaining observations, and hence with lower 
efficiency than for the pn camera. Nevertheless, we analyzed the mos data and found the 
resulting fluxes and parameters were quantitatively in agreement with those from pn data, 
given the current knowledge of cross- calibration uncertainties between the two instrumentso 
In this paper we report only the higher quality EPIC pn data. 

The data were analyzed with the XMM Science Analysis System (SAS) version T.l.clfl 
and the latest calibrations. Strong background flares can sometimes contaminate source 
events. To exclude possible flares, we extracted light curves from the entire field of view for 
events having energy > 10 keV. We then examined these light curves for flares. We defined 
bad time intervals to be when flares occurred, and excluded these intervals for all subsequent 
analyses. For all the XMM observations, we filtered a total of 20 ks of bad time intervals. 
Then we corrected the event times to the barycenter using the SAS barycen tool. 



2.2. RXTE observations 



We have observed AXP IE 2259-f586 regularly since 1997 with RXTE (see, e.g.. lGavriil fc Kaspi 
20021 ) . Our data were ob tained using the Proportional Counter Array (PCA) on board 



RXTE Jjahoda et al.ll2006h . The PCA consists of five identical and independent xenon/methane 
Proportional Counter Units (PCUs). We use our /?XTi? observations of IE 2259-1-586 to 
monitor its pulsed flux, and its frequency evolution using phase-coherent timing, and to look 
for bursts and pulse profile changes. 

For the purposes of this paper we analysed 193 observations that took place between 2001 
April 1 (MJD 52,000) and 2006 September 22 (MJD 54,000): 15 preoutburst observations, 
1 observation during the outburst, and 177 postoutburst observations. All 193 observations 



^See |http://xmm.esac.esa.int/docs/docuinents/CAL-TN-0018-2-6.pdf on the current calibration status 
of the EPIC cameras. 

^See http://xmni.esac. esa.int/sas/7. 1.0/ 
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with the exception of the two observations immediately following the outburst were taken in 
GoodXenonwithPropane or GoodXenon data modes. Both data modes record photon arrival 
times with 1 fis resolution and bin photon energies into one of 256 channels. To maximize the 
signal-to-noise ratio, we analysed only those events from the top xenon layer of each PCU. 
The remaining two observations were in event modes with a time resolution of ~125 /xs, a 
smaller number of energy channels, and no possibility of layer selection. For each of the 
observations we created barycentered light curves in the 2-10 keV band with 31.25 ms time 
resolution. 

We then folded each of the light cu rves using an ephemer is determined iteratively by 



maintaining phase coherence (see, e.g., iGavriil &: Kaspil |2002| ). We then used the folded 



pro files to calculate th e pulsed flux for each observation using both an rm s estimator (see. 



e^ 



2008 



Woods et al.ll2004l ) and an area estimator after baseline subtraction (see [Archibald et al. 
in preparation, for details). The results obtained using the two methods were consistent 
with each other. Here we only report the area pulsed flux because, while more sensitive to 
noise, it is the quantity of primary interest. 

To calculate the area pulsed flux for a given folded time series, we used the following: 

N 



PFarea ^ ^ {pi Pmin) / i 



where Pi refers to the count rate in the ith bin, is the number of phase bins, and Pmm is the 
average count rate in the off-pulse phase of the proflle, determined by cross- correlating with a 
high signal-to- noise ratio template, and calculated in the Fourier domain after truncating the 
Fourier series to six harmonics. Finally, we combined the pulsed flux numbers from each of 
two consecutive weeks into a single number, with the exception of the burst observation and 
the two observations that followed it, which remained unbinned. The results are presented 
in Figure [T]a. 



3. Analysis and Results 

3.1. Spectrum evolution 

Source spectra were extracted from circular regions of 32". 5 radius around the source 
center, using the barycentered, flltered event flle described in §2.1[ Background spectra were 
extracted from circular regions of 50" radius centered ~ 3' away from the source center. 
For the observations taken in Small Window Mode, we extracted single- and double-photon 
events and excluded events on or close to a bad pixel using the fllter expression "FLAG 
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= && PATTERN <= 4". In the Full Frame Mode observations, the source is highly 
off-center in the CCD image (Tabled]), and bad pixels were found close to the source center 
region. For these observations, the event list was filtered using the selection expression 
^XMMEA_EP to exclude only photons which fall directly on the bad pixels. However, we 
did not exclude photon events located adjacent to the bad pixel (which normally would 
be excluded by the expression FLAG = 0), because when there is a bad pixel close to 
the center of the source region, the effective area is evaluated more accurately with pixels 
around the bad ones taken into account by the SAS command arf gen [XMM help desk 2008, 
private communication). In order to avoid events that affected multiple pixels, we used only 
single events (PATTERN= 0) in the Full Frame Mode data. Event lists thus extracted 
were input to f tool grppha, which grouped the events by at least 25 photons per bin. A 
systematic uncertainty of 2% was also appended to the output spectra using grppha in order 
to characterize the current level of calibration accuracyjfl 



The resultant spectra were fitted in XSPEC 12.3.03 with the commonly used photoelec- 
trically absorbed blackbody plus power law model in the energy range 0.6-12 keV. Because 
the hydrogen column density Nh is not expected to be variable, we fixed this parameter for 
all the data sets and performed a joint fit. The goodness of fit is reasonable (see in Table 
[2]). The best-fit Nh is (1.012 ± 0.007) x lO^^cm'^. 

This value is consistent with that estimated from fittin g individual absorption edges o f 
elements O, Fe, Ne, Mg, and Si in the XMM RGS spectra (jPurant fc van Kerkwijkll2006l ). 
The other parameters were set free to vary and their best-fit values are presented in Table 
[21 The best-fit blackbody temperature, blackbody radius, and power law index are plotted 
versus time in Figure [H 

In order to look for correlations between spectral hardness and flux as observed i n 
other AXPs jRea et aLlboOsI : ICampana etaPboOTl : iTam et allboOSbl : IConzalez et aPboOsI ). 
we have looked for a correlation between hardness ratio and observed flux. We define the 
hardness ratio to be the ratio of 2-10 keV absorbed flux to 0.1-2 keV absorbed flux. We 
find the hardness ratio to be strongly correlated with the 2-10 keV absorbed flux (as shown 
in Fig. [2]a) in our observations. An anti-correlation between photon index and 2-10 keV 
unabsorbed flux is also seen, but has more scatter (as shown in Fig. This is likely 

because the photon index is not a perfect measure of spectral hardness, as it can be strongly 
influenced by the spectral flt at the low end of the band. 



''See http://xinm.esac.esa.int/docs/docuinents/CAL-TN-0018. 

Analysis) 

5 



[EPIC Status of Calibration and Data 



^See |http://heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
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3.2. Pulsed fractions 



We folded the 0.1-2 and 2-10 keV light curves of each XMM observation at the pulsar's 
period, deter mined usinK an e phemeris derived by phase coherent timing, from RXTE monitoring 
(Table [U see iDib et al.l 120071 for details). Each pulse profile was constructed by folding the 
photons into 32 phase bins. We measured area pulsed fiux of the XMM the same way we did 
for RXTE {see eqJT]), except t hat we used eight har monics instead of six when smoothing 
the light curves (for details see [Archibald et al.ll2008l in preparation). 



The measured area pulsed fractions are plotted in Figure [T]f. A possible correlation 
between the 0.1-2 keV area pulsed fraction and the 2-10 keV unabsorbed fiux is seen (Fig. 
[3l filled circles). A similar correlation was also found between the 0.1-2 keV area pulsed 
fraction and the 0.1-2 keV absorbed fiux. However, the correlation between 2-10 keV pulsed 
fraction and fiux is not significant (Fig. [3l open boxes). 

We also measured the rms pulsed fraction from the profiles to compare with the area 
pulsed fraction results. The 2-10 keV rms pulsed fractions are consistent with being constant, 
while the 0.1-2 keV rms pulsed fractions have some variance, but no significant trend or 
correlation with other parameters. The area and rms pulsed fractions are different by a 
factor that depends on the shape of the profile; as the pulse profile of IE 2259+586 did 
change temporari l y after the outburs t (from a simple double peaked profile to triple peaked; 



Kaspi et al 



2003 



Woods et al.ll2004j ). the different result is not surprising. 



3.3. Flux evolution 

We fit the unabsorbed fiuxes measured in our XMM observations after the outburst 
with a power law plus constant decay model, F{t) = Fb[(t — tg)/{l day)]" + Fg, where F{t) 
denotes the unabsorbed fiux, Ff, is the unabsorbed source fiux one day after th e onset of the 



outburst, Fq is the quiescent fiux and tg marks the glitch epoch MJD 52,443.13 (IWoods et al. 



20041 ). A good fit of x^i.^) = 0.66(5) (Fig. HI dashed line) was found. The best-fit power law 
index a = —0.69 ± 0.03 (Table [3]). The quiescent fiux level we found from this power law fit 
is (1.75±0.02) X 10~^^ ergss~^cm~^, considerably higher than that measured one week before 
the outburst [(1.59 ±0.01) x 10~" ergs s'^cm'^; Table[2]. We also fit the XMM unabsorbed 
fiuxes with an exponential decay plus quiescent level model, F{t) = Fpe"^*^*^^/'^ + Fg, where 
F(t) is unabsorbed fiux, Fp is the peak fiux, Fg is the quiescent fiux, r is the decay timescale 
and tg marks the glitch epoch. The fit is worse than that of the power law decay model 
but still acceptable, with xli^) of 1.08(5). The best-fit decay timescale r is 13.3 ± 0.7 days. 
Best-fit fiux decay parameters are presented in Table [31 
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We also fit power law and exponential models to the area pulsed flux of IE 2259+586 
measured by RXTE from 12 to 1649 days after the glitch. A power law model fits the 
data much better than the exponential model [xli^) = 1-18(69) for the power law model, 
xli^) = 1.65(69) for the exponential decay model; Table [3], which is evidence against the 
latter. The best-fit exponential decay timescale for RXTE data is 134 ± 15 days, an order 
of magnitude different from the ~ 13 day timescale found for the XMM data. 

The best-fit power law plus constant model for the evolution of the RXTE pulsed fluxes 
is different from that of the XMM total fluxes. This suggests that the 2-10 keV pulsed 
fractions were varying. In principle, we can check this with the pulsed fraction measurements 
we made with XMM (see § 13. 2p . Given the uncertainties on the XMM 2-10 keV pulsed 
fractions (Table [2]), as well as those of the best-fit evolution models (Table [3]), we find that 
the two are in agreement. 



Gotthelf fc HalpernI (120071 ) fit the spectrum of XTE J1810-197 using a double-blackbody 
model when studying that source's relaxation following its outburst. In order to com- 
pare the spectrum and evolution of IE 2259+586 to that of XTE J1810— 197, we also 
fit a photoelectrically absorbed double-blackbody model to IE 2259+586's spectra jointly. 
A double-blackbody model does not fit the spectra as well as the blackbody plus power 
law model (see Table [2] for details). The best-fit Nh for the double-blackbody model 
[(0.568 ± 0.03) X 10^^ cm~^] is smaller than that from our blackbody plus power law fit 
and is not consistent wit h the value measured model independently from RGS spectra 
[(1.12 ± 0.33) X 10^^ cm~^; iDurant &: van Kerkwiik|[2006l. but is consistent with the best-fit 



Nh [(0.5-0.7) X 10^^ cm-^] of CTB 109 measured bv ISasaki et all mm . 



Unabsorbed fiuxes obtained using the double-blackbody spectral model can also be 
fitted to a power law decay or an exponential decay model. The best-fit power law index is 
—0.73 ± 0.04, and the best-fit exponential timescale is 12.7 + 0.7 days (Table [3]), consistent 
with what we obtained using the blackbody plus power law spectral model. This indicates 
that our results for the decay parame ters are independent of the choice of spectral model. 



In the analysis of XTE J1810-197 bv iGotthelf fc HalpernI (120071 ) . they found that both of 



the two-blackbody components' fiux followed an exponential decay after XTE J1810— 197's 
2003 outburst. However, we find that the fiux of IE 2259+586's soft blackbody component 
measured from our fourth and fifth observations (only ~21 days after the outburst and glitch) 
were lower than that measured for the last five observations (see Table [2]for details). This fiux 
variation of the soft blackbody component therefore cannot be well fitted with an exponential 
or power law decay model. The temperatures of both the hotter and cooler components were 
also lower in the fourth and fifth observations than in the last five observations. The non- 
monotonic variation of the soft blackbody fiux and the two components' temperature are 
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different from what was observed by iGottfielf fc HalpernI (120071 ) and suggest that the double- 
blackbody model is not a reasonable representation of the spectrum of IE 2259+586. On 
the other hand, the spectral evolution from the blackbody plus power law spectral fit looks 
more reasonable. Using this spectral model, the blackbody radius in the first postoutburst 
observation was small compared to that of the preoutburst observations and was even smaller 
in the second and third postoutburst observations (Fig. [T]g?), suggesting that one or more hot 
spots formed after the outburst and were fading away in the next few months. In the last 
five observations, the blackbody radius was as large as the preoutburst value, suggesting that 
the putative hot spots had completely faded away and the thermal radiation then mostly 
came from the bulk surface of the neu tron star as i t did before outburst. Perhaps a more 
realistic spectral model such as that of iGiiver et al.l (120071 . l2008l ) could describe the spectral 
evolution of IE 2259+586 better, but such an analysis is outside the scope of this paper. 



Based on RXTE observations. IWoods et all (12004 ) found that the decay of IE 2259+586's 
2002 outburst consisted of two parts: a steeper power law decay in the first few hours, and 
a slower power law decay afterwards. They also found that the total energy released (2-10 
keV) in the slower decay was 2.1 x 10^^ ergs , which is much larg er than the total energy (2-60 
keV) released in the bursts (6 x lO^^ergs; iGavriil et al.ll2004l ). We also studied the slower 
decay, b y fitting a power law plus constant model, instead of the simple power law model 
used by IWoods et al.l (120041 ). The total rele ased energy, a ccord ing to our best-fit model, 
is roughly consistent with that calculated by IWoods et al.l (120041 ): we find 3 x 10^^ ergs 
(2-10 keV), assuming that the outburst will be over in 10000 days. However, based on our 
best-fit exponential model, the total energy released was somewhat smaller, ~ (3 — 4) x lO^'' 
ergs (2-10 keV). 



4. Discussion 



In this paper, we have presented a comprehensive study of the X-ray recovery of AXP 
IE 2259+586 following its 2002 outburst. Here we discuss the properties of this recovery, 
compare them with those of other magnetar outbursts, and consider how they constrain the 
magnetar model. 



4.1. Return to "Quiescence 



In our 2004 and 2005 XMM observations, the source's temperature and unabsorbed 
fiuxes were still higher than the preoutburst value (Fig. [T]). This suggests that the source 
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was not fully back to the preoutburst flux level. Our power law fit to the flux decay shows 
that the after-outburst quiescent flux level is (1.75 ± 0.02) x 10~^^ ergs, s~^cm~^, which is 
significantly higher than the preoutburst value [(1.59 ±0.01) x 10~^^ ergs s~^cm~^; Table [3]. 
Either the 2005 flux had still not returned to its quiescent level, or perhaps it had returned 
to quiescence but the flux just before the event was unusually low. Also possible is that this 
(and other) AXPs do not have well-defined constant quiescent fluxes, but have long-term 
flux variations. Indeed, there is evidence for some X-ray flux va riability in IE 2259+586 
over the years since its discovery in 1981 JSavkal fc Swanklll996h . Other AXPs also show 
variability on a variety of timescales (see iKaspil 120071 for a review) . 



4.2. Comparison with other Magnetar Recoveries 



It is useful to compare the behavior observed from IE 2259+586 with that of other 
magnetars. SGR 1900+14's f lux was found to fo llow a power law of index —0.713 ± 0.025 
after its 1998 August 27 flare JWoods et al.lboOlh FI This has been interpret ed as the cooling 



of th e magnetar outer crust following a sudden release of magnetic energy (ILyubarsky et al. 
2OO2I ). This model predicts a power law decay of index ~ —2/3. The flux of SGR 1627—41 
was found to decay following a power law of index ~ —0.47 since its 1998 source activation. 
Approximately 800 days after the source activation, SGR 1627— 41's flux suddenly declined 
by a factor of 10. This behavior is also we ll fitted by the crust cooling model, although with 
some fine tuning (IKouveliotou et al.l 120031 ). We fit the XMM 2-10 keV unabsorbed fluxes of 
IE 2259+586 with a power law plus constant model, and found the best-fit power law index 
to be —0.69 ± 0.03, close to that of SGR 1900+14, and that predicted by the model. This 
suggests that the IE 2259+586 outburst afterglow may also be explained by the diffusion of 
heat in the outer crust. 



The transient AXP XTE J1810— 197 exhibited an outburst in 2003. Ibrahim et al. 



(120041 ) found that the afterglow of the XTE J1810-197 outburst as observed by RXTE could 
be described by a power law decay model {F oc t~^) with P = 0.45 — 0.73. T his is similar to 



the behavior of IE 2259+586 and other SGRs. However, iGotthelf fc HalpernI (120071 ) found 
that, with more observations taken by Chandra from 2003 to 2006, the afterglow of the 
XTE J1810— 197 outburst actually followed an exponential decay of timescale 233.5 days. 
As we have shown in this paper, the pulsed and unabsorbed X-ray flux decay of IE 2259+586 



^Later the afterglow of the SGR 1900+14 Augus t 27 flare was fitted with a power law plus constant 
model ins tead of the single power law model used bv IWoods et al.l (|2001l) . and a decay index of ~ 0.9 was 
obtained (jFeroci et al.ll2003r ). 
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favors the power law decay model over the exponential decay. Perhaps the physical processes 
involved in the 2003 outburst of XTE J1810— 197 were different from those in 2002 outburst 
of IE 2259+586. 



4.3. Twisted Magnetosphere Model 



Thompson et al.l (120021 ) reported that, if there exists a global twist of the magnetosphere, 
the decay timescale r of this twist would be 



r 



4OA0^ 



pole \2/ -^NS \3 



1035ergss-i' 'lOi^G' 'lOkm' 



yr. 



(2) 



Woods et all fl2004 ) argued that, for IE 2259+586, the twist angle A0 should be ~ 10"^ 
rad. Thus, the predicted twist relaxation timescale of IE 2259+586 is several hours, which 
is coincident ly the timescale of the steeper flux decay observed at the beginning of the 
afterglow. 



However, iBeloborodov &: ThompsonI (120071 ) have shown more recently that this decay 
timescale is actually expected to be much larger than equation ([2]) suggests. This is because, 
in their model, the self-induction of the twisted portion of the magnetosphere accelerates 
particles from the stellar surface and initiates avalanches of pair creation which forms the 
corona. This corona persists in dynamic equilibrium, maintaining the electric current, as 
long as dissipation permits. The relevant timescale in this picture for the decay of a sudden 
twist is given by 



0.3( 



L 



X 



t)( 



(3) 



lO^^ergss ^^^GeV ^lOkm 
where Lx is the peak X-ray luminosity and e^e is the voltage along the twisted magneti c field 
lines and should nearly universally be ~1 GeV (see IBeloborodov fc Thompson! 120071 ). For 
IE 2259+586, we find r ~ 1.2 yr, given the peak luminosity Lx ~ 4 x 10^^((i/3 kpc) ergs s^^. 
Thus, the longer observed decay after t he initial steep decline may i ndeed correspond to the 
untwisting of a coronal flux tube in the IBeloborodov fc ThompsonI (120071 ) picture, although 



the predicted time scale is somewhat smaller than the observed time to return to quiescence. 
We note that the IBeloborodov fc ThompsonI (120071 ) model predicts a linear flux decline, 
in contrast to wh at we have observed for IE 2259+586 and what has been observed for 
XTE J1810-197 JCotthelf fc HalDernl l2007h. Moreover, in the ~5 yr of RXTE monitoring 
of IE 2259+586 prior to its 2002 outburst (iGavriil fc Kaspill2002l ). its pulsed X-ray luminosity 



in t he 2-10 keV band was roughly con stant at ~ 2 x 10^^ ergs s ^. This also is puzzling given 



the IBeloborodov &: ThompsonI (120071 ) prediction that if the time between large-scale events 
is longer than the decay time from the previous event, the magnetar should enter a quiescent 
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state in which the observed luminosity is dominated by the surface blackbody emission. Why 
should the "quiescent" blackbody emission from IE 2259+586 be a full order of magnitude 
larger than that from XTE J1810— 197, especially given the latter's much larger inferred 
magnetic field (1.7x10^^ versus 6x10^^ G)? This disparity in "quiescent," steady luminosities 
is even larger when considering AXP IE 1841—045, which has an apparently steady 2-10 keV 
luminosity of 1.4 x 10'^^ ergs s~^, and comparing with probable AXP AX 1845— 0258 , whic h 
has quiescent luminosity approximately 2 orders of magnitude smaller (ITam et al.l 120061 ) . 
Distance uncertainties may contribute but not on a scale that can significantly alleviate this 
problem. This remains an interesting puzzle in magnetar physics. 

The twisted magnetosphere or flux tube models generically predict that the flux and 
spectral hardness of magnetars in outburst should be roughly correlated due to increased 
scattering optical depth when the t wist is larger. How ever, a similar prediction for a 
fiux/hardness correlation was made by lOzel fc Giiveii (120071 ) in their thermally emitting mag- 
netar model, using a simple prescription for the magnetosphere and scatt ering geometry, with 
the latter stationary, i.e. invoking no variable magnetospheric twists. iGiiver et al.l (120071 ) 
found that their model could reproduce the existing data for XTE J1810— 197. We note 
that hardness-in t ensity correlations have n ow been observed for RXS J170849. 0—400910 
(ICampana et al.l 120071 ). IE 1048.1—5937 (ITam et al.l l2008bl ). and as we report, in our 
IE 22 59-1-586 XMM observations. It would be interesting to apply analysis of lOzel &: Giiver 
(120071 ) to these data, but it is outside the scope of this paper. 



4.4. Other Observed Recovery Properties 



The fact that the rms and area pulsed fractions remained largely constant while the 
blackbody radius (in the blackbody plus power law model) changed by a factor of ~2 (Fig. 
[1]) is worth considering, if the empirical blackbody plus power law spectrum model somehow 
resembles the real radiation mechanism. Pulsed fraction should generally decrease when the 
thermally radiating region on the star grows, provided that this region is not very small 
compared to the entire surface. Any realistic spectral model which takes radiative transfer 
in the atmosphere and scattering through the magnetosphere into account should be able to 
reproduce the observation in this regard as well. 

A clear anti-co rrelation between IE 1048.1— 5937 's pulsed fraction and u nabsorbed flux 
has been observed (iTiengo et al.ll2005l : iGavriil et al.ll2006l : iTam et al.ll2008bl ). However, we 
found no such correlation in the 2-10 keV band for IE 2259-1-586. On the contrary, its 0.1-2 
keV area pulsed fr actions seem to be correlat ed with both 0.1-2 and 2-10 keV unabsorbed 
fluxes (see Fig. E]). lOotthelf fc HalpernI J2007h found that XTE J1810-197 's pulsed fraction 
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measured between 2003 and 2006 after its outburst decreased with the decay of its flux, i.e. 
XTE J1810— 197's pulsed fraction is also correlated with flux. Thus, the striking anti- 
correlation between pulsed fraction and flux observed from IE 1048.1—5937 is clearly not 
universal. 

Finally, we note that the near-infrared flux decay of IE 2259+586 was foun d to follow 



a pow er law of index — 0.751;q 33 when fitted to a power law plus constant model ( Tam et al. 



2004 ). This decay index is close to what we found for the X-ray flux dec ay, thus confirming 
the reported correlatio n between near-IR and X-ray fiuxes postoutburstFI E^m et al. (2008a) 
and lWang et al.l (120081 ) showed that the nea r-IR fiux of IE 1048. 1—5937 do show correlation 
with X-rays at times of outbursts. However, Camilo et al.l (j2007l ) show that the near-IR fiux 
variation of XTE J1810— 197 is not simply correlated with X-ray fiux nor even monotonic 
postoutburst. Thus, the AXP picture with regard to near-IR variability is not yet fully clear. 
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Table 1. ZMM observation log for IE 2259+586. 



Name 


XMM 


Date 


Date 


pn on-times ^ Off-axis angle 


Frequencies ^ 




Obsid 


(MJD TDB) 


(YY-MM-DD) 


(fcsec) (arcmin) 


(s-i) 





0057540101 


52,296.791 


02- 


01- 


22 


8.5 


8.7 


0.1432871204(7) 




0038140101 


52,436.413 


02- 


06- 


11 


26.3 


2.0 


0.143287000(9) 


Obs3 


0155350301 


52,446.446 


02- 


06- 


21 


16.4 


2.0 


0.14328759(1) 


Obs4'i 


0057540201 


52,464.368 


02- 


07- 


09 


5.2 


10.7 


0.14328754(1) 


Obs5^ 


0057540301 


52,464.602 


02- 


07- 


09 


10.2 


10.3 


0.14328754(1) 


Obs6 


0203550301 


53,055.596 


04- 


02- 


20 


3.6 


1.9 


0.143286974(7) 


Obs7 


0203550601 


53,162.655 


04- 


06- 


06 


4.8 


2.0 


0.143286882(4) 


Obs8 


0203550401 


53,178.634 


04- 


06- 


22 


3.4 


2.0 


0.143286868(2) 


Obs9 


0203550501 


53,358.014 


04- 


12- 


19 


3.5 


2.0 


0.143286714(1) 


ObslO 


0203550701 


53,579.970 


05- 


07- 


28 


3.3 


1.9 


0.143286523(7) 



^On-times quoted reflect on-source times after filtering of background flares. 
''Frequencies are from contemporaneous RXTE observations. 



Observations taken before the outburst, MJD 52,443.13 (| Woods et all 120041 ). 

■^These three observations were taken in extended Full Frame Mode; all the others were taken in Small 
Window Mode. 
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Table 2. IE 2259+586's best- fit spectral parameters and pulsed fractions. 



Parameter'* ^ 


Obsl 


Obs2 


Obs3 


Obs4 


Obs5 


Obs6 


Obs7 


Obs8 


Obs9 


ObslO 








Blackbody plus 


power law model 










Nh (10^^ cm~2) 


1.012(7) 


1.012(7) 


1.012(7) 


1.012(7) 


1.012(7) 


1.012(7) 


1.012(7) 


1.012(7) 


1.012(7) 


1.012(7) 


kT (keV) 


0.37(1) 


0.406(2) 


0.510(4) 


0.48(2) 


0.49(1) 


0.400(7) 


0.400(5) 


0.400(7) 


0.405(7) 


0.400(7) 


r 


3.75(4) 


3.89(2) 


3.49(2) 


3.71(5) 


3.72(4) 


3.77(3) 


3.78(3) 


3.79(3) 


3.79(3) 


3.75(3) 


Fliix° 


1.15(2) 


1.29(1) 


3.45(3) 


1.95(5) 


2.03(4) 


1.51(2) 


1.48(2) 


1.48(2) 


1.47(2) 


1.45(2) 


Unabs Flux^ 


1.41(3) 


1.59(1) 


4.12(3) 


2.34(7) 


2.44(5) 


1.84(3) 


1.82(2) 


1.81(3) 


1.80(3) 


1.77(3) 


PL/BB ratio" 


1.8(3) 


1.2(2) 


1.6(2) 


2.6(4) 


2.1(3) 


1.6(2) 


1.5(2) 


1.7(2) 


1.5(2) 


1.7(2) 


Hardness ^ 


0.93(3) 


0.94(1) 


1.43(1) 


1.13(4) 


1.15(4) 


0.98(2) 


0.98(2) 


0.97(2) 


0.98(2) 


0.99(2) 




1.02(5800) 




(P = 


0.12)g 














Double blackbody model 


Nh (10^^ cm-^) 


0.568(3) 


0.568(3) 


0.568(3) 


0.568(3) 


0.568(3) 


0.568(3) 


0.568(3) 


0.568(3) 


0.568(3) 


0.568(3) 


Cooler kT (keV) 


0.362(5) 


0.372(2) 


0.390(3) 


0.330(7) 


0.335(5) 


0.371(3) 


0.380(3) 


0.371(3) 


0.370(4) 


0.371(3) 


Hotter kT (keV) 


0.77(4) 


0.82(1) 


0.86(1) 


0.74(3) 


0.73(2) 


0.85(3) 


0.94(2) 


0.89(3) 


0.85(3) 


0.86(3) 


Flux'= 


1.12(6) 


1.27(2) 


3.34(4) 


1.85(9) 


1.93(7) 


1.48(4) 


1.46(3) 


1.44(4) 


1.44(4) 


1.41(4) 


Unabs Flux<* 


1.26(6) 


1.42(2) 


3.69(4) 


2.05(11) 


2.15(8) 


1.65(5) 


1.63(4) 


1.62(5) 


1.61(5) 


1.58(5) 


HB/CB ratio^ 


0.8(1) 


0.6(1) 


1.2(2) 


2.1(3) 


2.0(3) 


0.7(1) 


0.53(7) 


0.64(9) 


0.7(1) 


0.7(1) 


Hardness ' 


0.91(5) 


0.92(1) 


1.37(2) 


1.07(6) 


1.09(4) 


0.96(3) 


0.95(2) 


0.94(3) 


0.95(3) 


0.96(3) 




1.11(5800) 




(P = 4.6 X 10-8)8 














Pulsed fractions 


PF(0.1-2 keV)' 


0.18(3) 


0.234(6) 


0.322(6) 


0.30(2) 


0.28(2) 


0.29(2) 


0.26(1) 


0.24(2) 


0.27(2) 


0.26(2) 


PF(2-10 keV)' 


0.23(5) 


0.30(1) 


0.339(9) 


0.33(4) 


0.36(3) 


0.30(3) 


0.33(2) 


0.34(3) 


0.31(2) 


0.29(3) 



"Numbers in parentheses indicate the Icr uncertainty in the least significant digit. Note that these uncertainties reflect the la error 

for a reduced of unity. 

''Best-fit parameters from a joint fit to all data sets. Nh in all data sets was set to be the same; other parameters were allowed to 
vary from observation to observation. 

"^(10—^^ ergs s^^cm^'^). Observed flux from both spectral components in the range 2—10 keV. 

''(10^"'^"' ergs s^^cm^^). Unabsorbed flux from both spectral components in the range 2-10 keV. 

'^The ratio of power law flux to blackbody flux in the 2-10 keV band (corrected for absorption). 

■''Spectral hardness defined as the ratio of 2-10 keV absorbed flux to 0.1-2 keV absorbed flux. 

"The probability for the xi to be higher than that was observed, assuming the model is correct. 

''The ratio of hot blackbody flux to cool blackbody flux in the 2-10 keV band (corrected for absorption). 

'The area pulsed fractions. 
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Table 3. Best-fit parameters for the IE 2259+586 flux decay 



Power law decay^ 


^,(10 ^-'^crgss ^cm ^) 


^5(10 -"^^crgss -'cm ^) 


a 




x'^/f 


XMM UF (BB+PL) 


(1.75 + 0.02) 


(5.40 + 0.21) 


-0.69 + 0.03 


3.31 


0.66 


XMMVF (BB+BB) ° 


(1.58 + 0.01) 


(5.07 + 0.22) 


-0.73 + 0.04 


1.22 


0.24 


RXTE PF d 


(0.14 + 0.06)cts S-IPCU-I 


(1.4 ± 0.1)cts s-lpCU-l 


-0.27 + 0.05 


81.3 


1.18 


Exponential decay" 


F,(10-"ergss-^cm-^) 


^^(lO^^crgss^icm-^) 


T (days) 




x'/^ 


XMMVF (BB+PL) 


(1.81 + 0.01) 


(2.97 + 0.06) 


13.3 ± 0.7 


5.39 


1.08 


XMM UF (BB+BB) <= 


(1.62 + 0.01) 


(2.69 + 0.05) 


12.7 + 0.7 


1.737 


0.35 


RXTE PF d 


(0.362 + 0.005)cts s-ipCU"! 


(0.46 + 0.03)cts S-IPCU-I 


134 + 15 


113.5 


1.65 



■^Power law decay model defined as F{t) = F(,{{t — tg)/(l day))" + Fq, where Fit) is unabsorbed flux, Fq is the quiescent 
flux, a is the power law index and tg is the glitch epoch MJD 52,443.13. 

''XMM unabsorbed flux decay measured using a blackbody plus power law spectral model. 

°XMM unabsorbed flux decay measured using a double-blax;kbody spectral model. 

''-RXTE area pulsed flux in units of cts s^ipCU^l. 

"Exponential decay model, defined as F{t) = Fpe~'*~*9)/'^ + Fq, where F{t) is unabsorbed flux, Fp is the peak flux, Fq is 
the quiescent flux, r is the decay timescale and tg is the glitch epoch MJD 52,443.13. 
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Fig. 1. — Spectral and pulsed fraction evolution of IE 2259+586 during and following its 
2002 outburst, (a) 2-10 keV area pulsed flux measured in RXTE monitoring observations; (fe) 
2-10 keV unabsorbed phase-averaged flux from XMM observations (all lower panels are also 
from XMM observations) ; (c) blackbody temperature (kT)^ (d) blackbody radius; (e) photon 
index; (/) 0.1-2 keV area pulsed fraction ( filled circl es) and 2-10 keV area pulsed fraction 
{open boxes). A distance of 3 kpc (IKothes et al.ll2002l ) is assumed to calc ulate the blackbod y 
radius. The vertical line denotes the 2002 ghtch epoch, MJD 52,443.13 JWoods et al.ll2004l ). 
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Fig. 2. — (a) Hardness vs. absorbed flux. Hardness ratio is defined as the ratio of 2-10 to 
0.1-2 keV absorbed flux. (6) Photon index vs. 2-10 keV unabsorbed flux. 
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Fig. 3. — The 2-10 keV area pulsed fraction {open boxes) vs. 2-10 keV unabsorbed flux ; 
0.1-2 keV area pulsed fraction (filled circles) vs. 2-10 keV unabsorbed flux. 
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Fig. 4. — Evolution of IE 2259+586's 2-10 keV unabsorbed phase-averaged flux 
(squares) measured from XMM observations and 2-10 keV area pulsed flux (cross) from 
i^XTE" observations following its 2002 outburst. XMM unabsorbed fluxes are in units of 
10~^^ergs s~^cm~^. RXTE area, pulsed fluxes are in units of count s~^PCU~^. The time 
axis is relative to the estimated glitch epoch (MJD 52,443.13). The solid line is the power 
law plus constant model fit to the RXTE area pulsed fluxes. The dashed line is a fit of the 
same model, although having different best-fit parameters, to the XMM fiuxes. See Table 
[3] for the best-fit parameters. The dotted line is the fiux level in 10~^^ergs s~^cm~^ ob- 
served with XMM owe week before the outburst. The uncertainty on this preoutburst fiux 
is approximately the width of the line. 



